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ABSTRACT 

Improvements on a statistical tropical cyclone (TC) track model in the western North Pacific Ocean are 
described. The goal of the model is to study the effect of El Nino-Southern Oscillation (ENSO) on East Asian 
TC landfall. The model is based on the International Best-Track Archive for Climate Stewardship ( IBTr ACS) 
database of TC observations for 1945-2007 and employs local regression of TC formation rates and track 
increments on the Nino-3.4 index and seasonally varying climate parameters. The main improvements are the 
inclusion of ENSO dependence in the track propagation and accounting for seasonality in both genesis and 
tracks. A comparison of simulations of the 1945-2007 period with observations concludes that the model 
updates improve the skill of this model in simulating TCs. Changes in TC genesis and tracks are analyzed 
separately and cumulatively in simulations of stationary extreme ENSO states. ENSO effects on regional 
(100-km scale) landfall are attributed to changes in genesis and tracks. The effect of ENSO on genesis is 
predominantly a shift in genesis location from the southeast in El Nino years to the northwest in La Nina 
years, resulting in higher landfall rates for the East Asian coast during La Nina. The effect of ENSO on track 
propagation varies seasonally and spatially. In the peak activity season (July-October), there are significant 
changes in mean tracks with ENSO. Landfall-rate changes from genesis- and track-ENSO effects in the 
Philippines cancel out, while coastal segments of Vietnam, China, the Korean Peninsula, and Japan show 
enhanced La Niha-year increases. 


1. Introduction 

There are many densely populated coastal regions that 
are susceptible to high fatalities and costly damage every 
year because of the frequent landfall of tropical cyclones 
(TC) that occur in the western North Pacific Ocean (WNP). 
Studies have shown that El Nino-Southern Oscillation 
(ENSO) strongly influences TC activity in the WNP. To 
better understand the TC landfall risk in this region, the 
relationship with ENSO must be considered. The most 
important change of WNP TCs with ENSO appears to be 
the shift in genesis location to the southeast during strong 
El Nino events and to the northwest during strong La 
Nina events (Chan 1985, 2000; Dong 1988; Lander 1994; 
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Saunders et al. 2000; Wang and Chan 2002; Chia and 
Ropelewski 2002; Camargo et al. 2007a). The shift is ex- 
plained by the location of the monsoon trough in the 
WNP shifting eastward during El Nino events (Lander 
1994) along with a reduction in vertical wind shear (Clark 
and Chu 2002), which both dictate favorable conditions 
for genesis. This shift results in changes in the tracks of 
the subsequent TCs: those that form in the southeast have 
longer lifetimes and tend to recurve northward as com- 
pared with those that form in the northwest that move 
straight westward onto land (Wang and Chan 2002). 
Eisner and Liu (2003) further studied the ENSO effect on 
tracks and implications for landfall, finding that strong El 
Niiio years correspond with northward recurving TC 
tracks, which are likely to make landfall on Japan, the 
Korean Peninsula, and northern China. During La Nina 
years, TC tracks move in a straighter westward path, 
putting the Philippines, southern China, and Vietnam at 
higher risk. The TC-track cluster analysis performed by 
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Camargo et al. (2007c) for the WNP demonstrates that 
the different genesis regions associated with El Nino (La 
Nina) years are also associated with longer (shorter) 
tracks. One important implication is that El Nino-year 
TC tracks in the WNP have more of an opportunity to 
reach higher intensities as they travel over the warm 
tropical ocean given a genesis location closer to the 
equator and farther eastward. El Nino-year TCs may also 
reach higher latitudes when they curve northward, before 
making landfall (Camargo and Sobel 2005; Chan and Liu 
2004; Camargo et al. 2007c). The last two effects are im- 
portant when making a landfall risk assessment in terms 
of landfall intensity and landfall location. 

The ENSO-track relationship has been primarily ex- 
plained by shifts in genesis location. One thing that has 
not been examined extensively is the direct effect of 
ENSO on tracks. For instance, independent of genesis 
location, will a TC tend to move in a different direction 
given different ENSO states? In this study, we in- 
vestigate this effect as well as the combined genesis- and 
track-ENSO effects to attribute the changes in landfall 
to each component. We fill a key gap in documenting the 
TC-ENSO relationship and its landfall implications. As 
much work has been done on improving ENSO pre- 
diction (e.g., Cane 1991; Fedorov et al. 2003; Goddard 
et al. 2010), a fuller understanding of the ENSO-TC 
relationship would allow for improved long-range prob- 
abilistic predictions of TC landfall activity (Vitart 2006; 
Vitart et al. 2007; Camargo et al. 2010). Here we present 
developments of a statistical track model that can 
simulate the way that East Asian landfall changes with 
ENSO. 

The advantage of using a statistical track model is the 
ability to translate what the TC activity-climate state 
relationships mean in terms of risk to humans. Obser- 
vational studies have used track shape (Eisner and Liu 
2003; Camargo et al. 2007b, c) or track density (Wang 
and Chan 2002; Wu et al. 2004; Zhang et al. 2012) changes 
as a proxy for risk of TC landfall because data become 
sparse when considering only the observed landfalls on 
small coastal regions. Data limitations increase when 
considering different climate states. Most other ENSO 
studies divide the observed data into subsets of less than 
10 years based on a threshold of an ENSO index. An 
important new aspect of our work is the use of stochastic 
simulations from a statistical track model that employs 
regression instead of creating subsets. Since the model is 
stochastic, its simulations are not identical to observed 
tracks. Our synthetic tracks are each unique but as an 
ensemble have similar statistical characteristics to the 
observed tracks. Another advantage of statistical models 
(Darling 1991; Drayton 2000; Vickery et al. 2000; James 
and Mason 2005; Emanuel et al. 2006; Rumpf et al. 2007; 


Hall and Jewson 2007; Graf et al. 2009; Yin et al. 2009; 
Yonekura and Hall 2011) is the ability to explicitly sim- 
ulate landfalls. Our track-modeling method allows for 
better landfall risk assessment for different ENSO states 
at a resolution of 100-km coastal segments instead of 
a country-scale assessment (e.g., Wu et al. 2004; Zhang 
et al. 2012). In effect, one could choose a specific segment 
or a city of interest instead of using a general risk level for 
a larger region or country. 

In previous work, Yonekura and Hall (2011, hereafter 
YH11) use a statistical track model to simulate the full life 
cycle of TCs in the WNP including an ENSO-dependent 
genesis component. Here, we describe the addition of 
ENSO dependence into the track-propagation component 
to the YH11 model. In YH11, the details of statistical 
model construction, validation, and simulation pro- 
cedures are presented. The genesis component depends 
on ENSO using local Poisson regression of the local 
genesis count time series on the July-October (JASO) 
Nino-3.4 index. The model was able to capture the geo- 
graphic characteristics of landfall in its simulations of the 
historical period, as well as replicate the spatial distri- 
butions of the observed track trajectories and track 
densities. 

One important aspect that was missing from the YH11 
model is the seasonality of WNP TCs. Chan (2000) has 
shown that the effect of ENSO on WNP TC tracks varies 
by month. For instance, in October of a strong El Niiio 
event, tracks will move more northward, avoiding Japan, 
because of midlevel westerlies around 40°N. In the same 
month for a strong La Nina, it is shown that tracks move 
westward in a straight path at lower latitudes because of 
easterly midlevel wind anomalies. Chia and Ropelewski 
(2002) documented the seasonality in the mean genesis 
location in the WNP, which moves from the southeast in 
March (5°N, 158°E) to the west in June (135°E), to the 
northeast in the boreal summer (17°N), and then back 
south-southeast in December (8°N). The monsoon 
trough moves in the peak typhoon season as well as with 
ENSO events, which explains the genesis-location shift. 
In the Camargo et al. (2007b, c) WNP track cluster 
analysis, it is also shown that different clusters have 
different seasonal occurrence rates as well as different 
seasonal periods of positive or negative correlation with 
ENSO. In this study, we include the effects of season- 
ality in the genesis and track components by using sea- 
sonally varying predictors, another improvement from 
the YH11 model. 

In the following section, we describe the methods used 
to incorporate ENSO and seasonal dependence into the 
model, run simulations, and discuss the model un- 
certainty. Then, in section 3, we evaluate the model, 
describe the effect of ENSO on TC tracks, and examine 
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simulations to attribute track and genesis effects on 
landfall rates. Section 4 concludes with a summary and 
discussion of the results. 

2. Methods 

a. Data 

As in YH11, the TC track data used to construct the 
statistical track model are from the International Best- 
Track Archive for Climate Stewardship (IBTrACS; 
Knapp et al. 2010) database of 6-hourly storm positions. 
The storms from 1945 to 2007 are used for consistency 
with YH11. 

To define the state of ENSO for a given year, the 
JASO Nino-3.4 index (anomalies from the 1950-79 
monthly climatology; Barnston et al. 1997) is used and 
will be hereafter referred to as the ENSO index. We 
acknowledge the recent studies of the central equatorial 
Pacific warming type of ENSO — an El Nino Modoki 
event (Ashok et al. 2007) — and would like to distinguish 
here that we are considering the conventional El Nino 
phenomenon defined by the first EOF from equatorial 
Pacific sea surface temperature (SST). The Nino-3. 4 and 
Nino-3 indices are both well correlated with the princi- 
pal component associated with the first EOF, and there 
is little difference in model results when using either 
index. 

Seasonally varying predictors are used in both the gen- 
esis and track-propagation component of the model. The 
main purpose of adding SST as a parameter in the genesis 
component and midlevel winds in the track-propagation 
component is to include the effect of seasonality on WNP 
TCs, thus we use climatological predictors. The interaction 
term between the climatological predictor and ENSO in- 
dex represents any existing correlation between the two 
that affects the TC model component, for example, if 
ENSO state affects seasonality or if the season affects the 
ENSO impact on TC activity. ENSO has known effects 
on the atmospheric circulation and SST in the WNP. At 
the midlevel troposphere, an El Nino event is associ- 
ated with a cyclonic anomaly over Southeast Asia, which 
enhances the westerly jet along the coast (Wang et al. 
2000). WNP SST, especially in 5°-20°N, 130°-155°E, is 
negatively correlated with the Nino-3.4 index though its 
variability is considerably less than the index region 
(Chan and Liu 2004). We believe that our model captures 
these effects on TCs indirectly through the ENSO index 
predictor. We assume the physical environmental changes 
are implicit in the observed TC behavior. 

For the genesis component, we use the monthly annual- 
cycle climatology of SST from the National Oceanic 
and Atmospheric Administration (NOAA) extended 


reconstructed SST, version 3 (ERSST V3), data (Smith 
and Reynolds 2004) provided by the NOAA/Office of 
Oceanic and Atmospheric Research/Earth System Re- 
search Laboratory/Physical Sciences Division (OAR/ 
ESRL PSD) in Boulder, Colorado (the data can be found 
at http://www.esrl.noaa.gov/psd). In the track-propagation 
component of the model, the 500-hPa daily annual- 
cycle climatology of zonal winds is used from the Na- 
tional Centers for Environmental Prediction-National 
Center for Atmospheric Research (NCEP-NCAR) 
reanalysis for the years from 1948 to 2007 (Kalnay et al. 
1996). 

The choice of 500-hPa zonal winds as a predictor for 
the track component is motivated by Chan and Gray 
(1982), who studied which wind levels most influence 
track displacement and speed. They found that the 
midlevel winds (700, 600, and 500 hPa) have the most 
influence. Chan (2000) then used 500-hPa wind anom- 
alies to explain changes in WNP TC tracks associated 
with different ENSO phases. Wu et al. (2004) made 
composites of the different ENSO-state 500-hPa winds 
to show that changes in landfall for El Nino and La Nina 
relate to the 500-hPa subtropical ridge. More recently, 
Zhang et al. (2012) looked at the 500-hPa wind-field 
anomalies and found mechanisms for ENSO-TC track 
effects. Looking at the NCEP-NCAR reanalysis data, 
the 500-hPa winds in the WNP follow a seasonal cycle. 
In December-February (DJF), there are strong west- 
erlies in the midlatitudes poleward of 20°N and weaker 
easterlies in the tropics. In March-May (MAM), the 
band of westerlies moves northward, especially in the 
eastern part of the basin. Then, in June-August (JJA), 
the equatorial easterlies strengthen, and there is also 
more northward motion toward the East Asian coast. 
By October, the midlatitude westerlies begin to shift 
back southward to 20°N. Here, we choose to use the 
500-hPa winds as a first effort to incorporate a physics- 
based parameter in track-propagation modeling. 

b. Genesis 

The genesis component of the model simulates the 
number of TCs per year and the day of year and for- 
mation location of each TC, all of which may be sensitive 
to ENSO. The total number is simulated by a random 
draw from a WNP-wide Poisson distribution whose rate 
is determined by Poisson regression of WNP-wide an- 
nual counts on the ENSO index (Sabbatelli and Mann 
2007; YH11). The determination of the location of 
a genesis event is as in YH11: local Poisson regression is 
performed using annual formation count data inside 
data circles whose size is determined by out-of-sample 
likelihood maximization. The spatial field of Poisson 
rates is normalized to create a 1° X 1° spatial probability 
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Fig. 1. Local Poisson rates calculated from the genesis model for El Nino (EN) and La Nina (LN) years (ENSO index values of +2 and —2, 
respectively). Select months that specify a monthly climatological SST are shown: (a) June EN, (b) August EN, (c) October EN, (d) June LN, 
(e) August LN, and (f) October LN. Red contours outline regions with highest expected genesis rates. 


density function (PDF), and a random draw from the 
PDF determines the location for each simulated TC. 
Further details can be found in YHlf . 

We have added seasonality to the genesis model of 
YFI11. We include the local monthly annual-cycle cli- 
matology of SST as an independent predictor variable in 
the regression. We randomly select a day of year, and 
the SST value of that day fixes a Poisson rate for 
a formation event. Then we sample the probability and 
repeat until an event is realized. An ENSO-SST interaction 
term is also used as a predictor to account for the effect 
of yearly ENSO state on seasonality. Equation (1) 
shows how to calculate a local Poisson rate A: 

A (r,a,m) = exp[/3 Q (r) + )S ENSO (r) X ENSO(fl) 

+ /3 SST (r) X SST(r,/n) + /^ensoxsstW 

X ENSO(fl) X SST(r, 777)] , ( 1 ) 

where r stands for a grid location, a specifies a year, 
777 specifies a month, each [3 is a local regression coeffi- 
cient at 7-, ENSO(fl) is the ENSO index for year a, and 
SST(t\ 777) is the SST at location r and month m. As in 
YH11, the shift in dominant genesis region to the south- 
eastern part of the WNP for El Nino years and to the 
northwest for La Nina years is observed for different in- 
puts of ENSO(a). shown in Fig. 1. The shift is most 
prominent between June and October, when 76% of 
WNP TCs occur. In Fig. 1, the genesis Poisson rate fields 
for June, August, and October in different ENSO states 


are shown. The subtle seasonal shifts of high-Poisson-rate 
regions in Fig. 1 agree with the seasonal cycle of mean 
genesis location from Chia and Ropelewski (2002). Our 
model now captures both ENSO and seasonal effects on 
WNP TC genesis. 

c. Tracks 

The track component of the model simulates the 
6-hourly increment of a TC position as a function of 
location, day of year, and ENSO state. The basic 
methodology is identical to that described in YH11, but 
here we have added seasonality and ENSO dependence. 
Our incorporation of ENSO and seasonality is guided 
by the work of Chan and Gray (1982) and Chan (2000), 
who analyzed ENSO-track relationships month by 
month and attributed the ENSO effect to anomalies in 
the 500-hPa winds. They demonstrated that the effect 
of ENSO on TCs in the WNP evolves through the season. 
Thus, the ENSO index and an annual cycle of 500-hPa 
winds are used as parameters to incorporate seasonality 
in track propagation. The wind data come from daily 
climatological zonal winds interpolated to a 1° grid. We 
found the meridional winds to be insignificant as a pre- 
dictor based on a two-tailed Student's t test. The merid- 
ional winds are not used to avoid modeling a relationship 
with a high degree of uncertainty, which could introduce 
a model bias. The lack of significance may be due to the 
meridional displacement being governed more by TC size 
and latitude. 
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Fig. 2. Local regression coefficients used in the track-propagation model component in the zonal ( x ) and me- 
ridional (y) directions. White areas indicate where the regression coefficient is not significant, (a) /Ienso* and 
(d) /Jensov are the coefficients that relate the mean track to the current-year ENSO index, (b) fi u soox relates the zonal 
mean track component to the daily climatological 500-hPa zonal winds, (c) /3ENSox u 5ou.r accounts for the correlated 
components of the ENSO index and 500-hPa zonal winds that effect the zonal mean track displacement. 


The IBTrACS 6-hourly zonal and meridional dis- 
placement values are linearly regressed locally on their 
respective predictors. The distance weight determining 
the locality of regression is a Gaussian kernel with 
a length scale of 300 km. The length scale is determined 
objectively by out-of-sample likelihood maximization. 
At each 1° grid point, we archive values of the regression 
coefficients from the equations: 

dX(r, a, d) = P 0x (r) + p ENSO (r) X ENSO(a) 

+ ^u500.v( r ) X M 50()( r ’ + ^ENSOXhSOOx^) 

X ENSO(a) X u 500 (r, d ) and 

(2) 

dY(r,a) = p 0y (r) + P EmQy (r) X ENSO(a), (3) 

where dX(r,a,d ) and dY(r,a) are the mean displace- 
ments in the zonal and meridional directions, respectively, 


at grid point r, year a , and day d. ENSO(a) is the ENSO 
index for the year specified by year a, and i/soo( r < d) is 
the 500-hPa zonal wind at grid point r and day d from 
a daily climatology. Note that ENSO is only included 
as an effect for the mean track increment. In full sim- 
ulations TCs propagate according to Eqs. (2) and (3) in 
the mean, and the residuals about the mean are mod- 
eled stochastically, as described in YH11: a bivariate 
anisotropic correlated normal, with lag-1 autocorrelation 
for the standardized errors. 

The resulting fields of coefficients from regression 
Eqs. (2) and (3) are shown in Fig. 2, in which the con- 
tours are color filled where the coefficients are signifi- 
cantly different than zero at the 95% confidence level 
according to a two-tailed Student’s t test at every grid 
point. ENSO significantly affects TC mean tracks. In 
Fig. 2a, the Penso.x regions of negative values near 20°N, 
150°E and 45°N, 150°E indicate that positive (negative) 
ENSO tends to enhance (inhibit) the westward motion 
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of TCs in an important region for landfall consideration. 
In this region, as well as in the rest of the basin, positive 
(negative) ENSO tends to inhibit (enhance) the north- 
ward motion of TCs (Fig. 2d). Figure 2b shows 
which not surprisingly shows a positive relationship 
between the midtropospheric zonal steering wind and 
TC motion. Finally, /3 EN sox„50(u (Fig. 2c) reveals 
some coupling between ENSO and 11500 , particularly 
at midlatitudes. For positive (negative) ENSO the TC 
motion enhances (inhibits) the relationship between 
TC motion and the seasonal steering winds. 

d. Lysis 

Lysis, or death, of each TC is determined by com- 
paring a random uniform draw to a predetermined PDF 
of historical lysis probability at every track step. Land 
and ocean lysis is strictly separated according to a 0.25° 
land-ocean mask. See YH11 for details. 

e. Simulations 

Our goal is to estimate the effects of ENSO on East 
Asian TC landfall via genesis and track propagation. 
Separation of the two ENSO-landfall effects is in- 
vestigated by running 10 000-yr simulations at stationary 
values of the ENSO index using different model ver- 
sions. Here, we present the simulations using the ENSO 
index values of +2 and —2 anomaly values from the 
1950-79 base period. The 10 000-yr duration is required 
to converge on mean landfall rates at high geographic 
resolution. The four model versions are 1) no dependence 
on ENSO; 2) genesis dependence on ENSO (“Genesis- 
Only”); 3) track dependence on ENSO (“Track-Only”); 
and 4) genesis and track dependence on ENSO (“Both”). 

We also evaluate each of the four model versions 
against historical landfalls. For this purpose, a separate 
set of simulations is performed: 500 realizations of the 
historical period 1945-2007 for each model version using 
the historical time series of the ENSO index. The mean 
and spread of each ensemble is compared to the his- 
torical landfall rates. The model is unbiased to the extent 
that the historical landfall rates appear as a typical 
member of the 500-member synthetic set. 

f. Uncertainty and significance 

There are two types of uncertainty to document. One 
is the uncertainty due to the stochastic nature of our 
simulations. For example, a single 1-yr simulation in a + 2 
ENSO state will have landfalls along the coast that differ 
from those of another 1-yr +2 ENSO simulation, 
mimicking the chaotic nature of TCs. When we per- 
form a multiyear simulation, we compute the mean 
value of a diagnostic (e.g., landfall rate) and the spread 
about the mean. How well do we know the mean? We 


perform long simulations for a fixed ENSO state to 
reduce the uncertainty of the mean due to finite simu- 
lation length. We find that 10 000-yr simulations are 
required to converge on mean model landfall rates at 
high resolution along the coast. We cannot, however, 
reduce the uncertainty due to the finite historical data 
on which the model is constructed. This is the second 
type of uncertainty, and it is crucial to document this 
uncertainty to say whether the relationships we un- 
cover between ENSO and a mean diagnostic value are 
significant. 

To determine the second type of uncertainty, we 
perform a generalized jackknife significance test (Tukey 
1958). We reconstruct the genesis and track models 100 
times, each time using a random —80% of the IBTrACS 
years (50 years out of 63). For each of the 100 jackknife 
members, we run 10 000-yr simulations at extreme 
ENSO states and compute the mean of the diagnostic 
(e.g., landfall rates along the coast). We then take the 
difference between +2 and —2 ENSO mean values for 
each jackknife member. If at least 97.5% of the dif- 
ference values are all positive or all negative, then we 
call the ENSO effect significant for that diagnostic. 

3. Results 

In this section we first evaluate the model by com- 
paring historical simulations to direct observations. We 
then examine the impact of ENSO on WNP TC tracks 
and landfall, first by illustrating the ENSO effect on 
mean tracks, and then by performing full stochastic TC 
simulations at stationary strong ENSO states and ex- 
amining the landfall rates. Landfall rates from three 
different model versions are used to isolate the effects of 
ENSO on landfall via genesis and tracks. 

a. Model evaluation and comparison 

We now evaluate the model by comparing simulation 
results to observations using two diagnostics: 1) life 
spans at extreme ENSO states and 2) landfall rates over 
the historical period. We calculate TC life spans from 
IBTrACS and find an average of 9.3 days for years when 
the ENSO index is greater than 1 (El Nino) and 6.8 days 
when the ENSO index is less than —1 (La Nina). Using 
the “Both” model version (ENSO-dependent tracks 
and genesis), we obtain a mean life span of 8.6 days 
during El Nino and 6.8 days during La Nina, in close 
agreement with the direct observations. Repeating the 
calculation with Genesis-Only and Track-Only model 
versions indicates that the ENSO influence on genesis 
causes the change in life span. Genesis-Only has 
a mean life span of 8.3 (6.3) days for El Nino (La Nina). 
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DISTANCE ALONG COAST IN 1000KM 


Fig. 3. The ensemble-mean landfall rates for each 100-km segment from simulations of the 
historical period (1945-2007) are shown with the solid red line. The red shading indicates the 
inner 95% spread of landfall rates over the 500-member ensemble. The Both model configu- 
ration was used, having ENSO-dependent genesis and tracks. Landfall rates from IBTrACS are 
in black. Letters highlight where the IBTrACS landfall rates fall outside the model inner 95% 
spread, corresponding to the letters in the inset. The map inset shows the coast defined by 
100-km segments used in the landfall diagnostic, separated into larger regions: main coast, 
Japan, and the Philippines. The dashed lines correspond to the dashed lines in the main panel. 


Track-Only has a mean life span of 7.5 (7.7) days. La 
Nina TCs originate, on average, closer to the coast and, 
thus, have a shorter life span before making landfall. 
We note that others (Wang and Chan 2002; Camargo 
and Sobel 2005) find the same sign in ENSO life-span 
change, though they use different ENSO thresholds. 
Camargo et al. (2007c) calculate the median lifetime 
for track types, or clusters, that occur more during El 
Nino states. Their result is higher than our average El 
Nino lifetime by about 2-4 days. The median lifetime 
for their cluster occurring most with La Niiia has about 
the same value as our La Niiia average. 

To evaluate the landfall characteristics of the model, 
we perform 500 full stochastic simulations of the his- 
torical period 1945-2007. We define a landfall as when 
a track crosses a 100-km coastal segment from ocean to 
land for the East Asian regions shown in the Fig. 3 inset. 
Figure 3 shows the average annual landfall rates per 
100 km over the 500-member ensemble, the inner 95% 
spread about the mean, and the rates directly computed 
from IBTrACS. Landfall rates vary enormously by re- 
gion, depending on the orientation of the coast with 
respect to mean tracks and the proximity of the coast to 
regions of high activity. Peak landfall rates occur in 
Southeast Asia, southern China, and the Philippines. 

The model is in close agreement with observed landfall 
rates. In most regions the historical landfall rate falls in- 
side the 95% spread about the mean, indicating no model 
bias. There are, however, a few scattered regions of bias 
that are identified with letters that correspond to the 
coastal segments highlighted on the Fig. 3 map inset. Just 
north (Fig. 3a) and south (Fig. 3b) of the Vietnam-China 
border, there are local landfall peaks where the model is 
biased low. Small regions of low bias also occur on the 


southeastern coast of Japan (Fig. 3c) and the western 
coast of the northern Philippines (Fig. 3d). The eastern 
coast of the southern Philippines (Fig. 3e) displays a small 
region of high bias. Despite these small biases, overall 
the analysis confirms that the model is able to capture 
the geographically complex behavior of TC landfall in the 
WNP. Furthermore, the inclusion of ENSO reduces 
model bias compared to a model version with no ENSO 
sensitivity (not shown). One possible reason for the bias is 
that the lysis model does not depend on the intensity of 
the storm. When an intensity component is added to the 
model, lysis will be determined more realistically by con- 
sidering intensity. Another possible explanation is that the 
landfall rates from 500 simulations of the 63-yr historical 
period may not be sufficient for the model to simulate the 
full range of possible landfalls. Flowever, it could also be 
that the “biases” are spurious. In any 63-yr realization of 
the stochastic process, one expects several local landfall- 
rate excursions beyond the 95% band about the mean, if 
enough regions are queried. 

b. ENSO effect on mean tracks 

The modulation of TC tracks in the WNP due to 
ENSO is examined by looking at mean tracks for ENSO 
index values of +2 and —2. Mean tracks are the trajec- 
tory followed as a result of the local regression of Eqs. (2) 
and (3), without any stochastic modeling of residuals. 
This is, in effect, the average trajectory that a TC would 
follow as it moves through the basin starting from 
a specified location and date. We compare mean tracks 
that originate from the same point to exclude the in- 
direct influence of ENSO on tracks via shifts in genesis 
site. 
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August 15 



100°E 140°E 180°E 


October 15 







Fig. 4. Local mean-track displacement vectors for ENSO-index values of +2 (red) and —2 (blue) shown for 
(a) 15 Aug and (b) 15 Oct, which each have distinct climatological daily 500-hPa zonal winds. Vectors are only 
shown where the local regression coefficients are significant. (c),(d) 500-hPa wind vectors from the daily climatology for 
each respective day. (e),(f) The difference between red and blue vectors from (a),(b), respectively, which highlights the 
regions where the local mean displacement changes with ENSO state. The units for the vectors in (a),(b),(e), and (f) are 
kilometers per step, where “step” refers to the 6-h steps of the model. 
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Fig. 5. Mean tracks initiated on (a)-(c) 15 Aug and (d)-(f) 15 Oct from numbered genesis locations with ENSO 
indices of +2 (red lines) and —2 (blue lines) for extreme El Nino and La Nina years, respectively. The dashed 
portions of the red and blue lines indicate where the 10-day mean tracks were extended to reach the diagnostic 
line. The dashed black lines are used for a diagnostic that measures the difference in crossing location between 
ENSO states. The asterisks by genesis point numbers indicate that the resulting mean tracks have a significant 
shift, and the red (blue) number colors indicate that the genesis site is more likely to occur during an El Nino 
(La Nina) year. 


Figures 4a and 4b show the local mean displacement 
vectors lor opposite ENSO states at two dates. The 
vector differences are highlighted in Figs. 4e and 4f and 
can be interpreted as showing the tendency of mean 
tracks in El Nino years relative to La Nina years. There 
are several key differences between mean-track vectors 
for the two ENSO states that are relevant to TC landfall 
risk, which are present in both months. At low latitudes 
(below 15°N), El Nino (red) vectors have weaker west- 
ward components than La Nina (blue) vectors, resulting 


in a slower shuttling of TCs toward the Philippines and 
Vietnam. South of Japan, red vectors exhibit a stronger 
northward component. West of Japan and around the 
Korean Peninsula, blue vectors have stronger northeast- 
ward magnitudes. Examining the 500-hPa wind from the 
daily climatology for the two dates in Figs. 4c and 4d 
shows the expected steering winds with easterlies in the 
tropics and stronger westerlies in the midlatitudes. 
Differences between August and October winds here 
include a stronger northward component around the 
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Table 1. Mean-track significance testing results for each genesis 
point and starting date. An “M” (“F”) signifies that the mean track 
shift from opposite ENSO states was significant for a starting date 
on the fifteenth (first) day of the given month. Asterisks indicate 
months when the ENSO effect on mean tracks is of the opposite 
sign from the other months. 




Genesis point 



1 2 3 

4 5 6 7 8 9 

10 11 

July 

M, F F 


M 

August 

September 

October 

F F 

p * * 

* * * 

M M, F 

F 

M* * F 

M 


Philippines in August, stronger easterlies in the tropics in 
August, and strong westerlies reaching lower latitudes in 
October; all differences are reflected in the mean track 
behavior in Figs. 4a and 4b. 

To better visualize the cumulative effect of these 
space- and time-dependent ENSO dependencies, we 
launch mean tracks of 10-day durations in +2 and —2 
ENSO states from a set of fixed locations in July, August, 
September, and October. Figure 5 shows the tracks for 
August and October, while Table 1 additionally docu- 
ments results for July and September. As a diagnostic of 
the ENSO effect we use the location where tracks cross 
lines of constant latitude and longitude. The ENSO 
difference in the location (degrees of latitude or lon- 
gitude) is deemed to be significant or not according to 
our jackknife uncertainty tests described above. An 
asterisk by the genesis point in each figure indicates 
that the ENSO effect is significant at the 90% confi- 
dence level. Figure 6 provides examples of the jack- 
knife test for the case of an insignificant ENSO effect 
and a significant ENSO effect. 

These 10-day tracks can be divided into two broad 
categories, those that recurve, ultimately heading north- 
eastward, and those that remain on a comparably straight 
northwestward trajectory. We find that both track cate- 
gories occur in both ENSO states, and that there is little 
variation from one category to another through the sea- 
son from a given genesis site. Genesis points 1-5 had 
straight-moving tracks through the whole season, while 
genesis points 6-11 had recurving tracks through much or 
all of the season. The only variations occurred for genesis 
points 6 and 10, which had straight -moving tracks in July 
and recurving tracks in August-October. 

ENSO influences the trajectory of both the straight- 
moving and recurving tracks. We summarize this in- 
fluence by calculating the shift in location that mean 
tracks cross specified lines of longitude or latitude 
(dashed black lines in Fig. 5). The sign of mean track 
shifts with ENSO changes through the season for the 
straight movers but stays constant for recurving tracks 




Fig. 6. Illustration of the significance test for the mean track shift 
between two ENSO states across diagnostic lines from Fig. 5. (a) 
The distribution of mean track shifts from a genesis site that results 
in an insignificant track shift: the inner 90% spread (shaded) 
bounds zero, (b) A similar distribution from another genesis site 
that produces a significant track shift. 

(Table 1, asterisks). The ENSO shift for straight-moving 
tracks from genesis points 1-5 (Fig. 5) reverses from 
August to October: in August, El Nino (La Nina) mean 
tracks move more westward (eastward) and northward 
(southward), and the effect reverses in October. Most of 
the shifts of either sign are not significant, though for 
select genesis dates there is a significant ENSO shift for 
points 1, 2, 4, and 5 (Table 1). Most of the recurving 
mean tracks from genesis points 6-11 propagate far- 
ther westward during El Nino before recurving. How- 
ever, only tracks from points 7 and 11 have significant 
shifts. 

The straight-moving tracks pose the biggest landfall 
threat for the Philippines, Vietnam, and southern China. 
The shifts with ENSO result in changes in landfall lo- 
cation within these regions. However, the ENSO effect 
reverses throughout the season resulting in some can- 
cellation of the net annual ENSO-track influence for 
these tracks. Recurving tracks threaten the coasts of 
Japan, the Korean Peninsula, and northern China. The 
El Nino westward shift of recurvaturc point enhances 
the chance of a landfall. 
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Before running full stochastic simulations, we first 
consider the combined landfall effect of ENSO on 
genesis and tracks by using the mean tracks and the 
known shift in genesis region. On Fig. 5, genesis points 
1, 2, 3, 7, and 11 have a higher likelihood of genesis during 
an El Nino year (red numbers), while genesis points 4, 5, 8, 
and 10 have a higher likelihood during a La Nina year 
(blue numbers) determined from Fig. 11 in YFI11. 
Looking at combined effects by following the red (blue) 
tracks for likely El Nino (La Nina) genesis points, results 
are similar to those of Eisner and Liu (2003): landfall rates 
increase in the northern Philippines and Japan during an 
El Nino and increase on Vietnam, China, the Korean 
Peninsula, and the Philippines during a La Nina. 

c. ENSO and landfall rates 

We now turn to landfall rates using the fully stochastic 
model. The results for a neutral state (i.e., ENSO index 
of 0) are not shown, but the neutral landfall rates fall in 
between the two ENSO states shown for most of the 
East Asian coast. As shown by YH11, the Genesis-Only 
model produces higher landfall rates everywhere 
along the coast during La Nina (Fig. 7a). The genesis- 
ENSO effect on landfall rates is most pronounced on 
the main coast (defined in the Fig. 7c inset) near the 
China-Vietnam border, where there are higher land- 
fall rates during La Nina (blue) than during El Nino 
(red). The result is expected: the region of primary 
genesis is shifted to the northwestern WNP during La 
Nina, closer to the coast. The genesis-ENSO effect is 
less pronounced for regions of low landfall rates, such 
as Indonesia-Malaysia and Kamchatka (not shown). 

The Track-Only model shows more spatially varied 
ENSO effects at landfall (Fig. 7b). The ENSO-track 
effect is similar to the ENSO-genesis effect for southern 
China, Vietnam, and the northern Philippines, having 
higher rates during La Nina. However, Japan, the 
southern Philippines, and Thailand have higher rates 
during El Nino, with the most pronounced increase in 
the Philippines. There is no longer a difference in the 
Korean Peninsula and northern China. Furthermore, 
landfall rates overall are higher using Track-Only than 
Genesis-Only. 

The landfall rates resulting from the inclusion of both 
track- and genesis-ENSO effects (the Both model) 
closely resemble those of the Genesis-Only model. 
Figure 7c shows that the ENSO effect on China, Vietnam, 
and the northern Philippines is enhanced because of 
the genesis and track effects combining. There is a less 
pronounced effect on Japan and little effect on the low- 
landfall regions. For the southern Philippines, there is 
virtually no difference between ENSO states, showing 
that the track effect offsets the genesis effect for this 


area. The red and blue shading show the 95% confi- 
dence, or uncertainty, on the mean landfall rates from 
the generalized jackknife significance test. In many 
regions, the jackknife uncertainty overlaps for the two 
ENSO states. 

To further highlight the significance of the ENSO ef- 
fect for the Both model, the difference in landfall rates 
from El Nino to La Nina is plotted in Fig. 7d, including 
the jackknife uncertainty of this difference. The regions 
where the uncertainty excludes zero have a significant 
ENSO effect at the 95% confidence level. The western 
and southern coasts of the Korean Peninsula, parts of 
China, part of Vietnam, southern Japan, and small points 
on the northern tips of Japan and the Philippines meet 
this condition (highlighted in blue on the Fig. 7c map 
inset). In these regions we can say with 95% confidence 
that there are more landfalls on average during La Nina 
than during El Nino. These results differ slightly from 
Saunders et al. (2000), who use a different type of sta- 
tistical model to determine the ENSO effect on East 
Asian landfall. They find that the only significant differ- 
ence is a La Nina increase in landfall in the Philippines 
(around Manila) and Vietnam. The former contrasts with 
our result, though it is in an area where our model ex- 
hibits a low bias. Saunders et al. (2000) also find an in- 
crease in landfall on Japan for El Nino and an increase in 
landfall on South Korea for La Nina, although it is not 
a significant change in their study. Our results agree 
well with those of Ho et al. (2005), where the frequency 
of TC passages in different observed ENSO years is 
examined. Their study did not find the effect to be as 
significant, possibly because they examine only 7 years 
of each ENSO state based on a threshold of ENSO 
index. 

Track density is now examined to help interpret these 
landfall-ENSO results. Density is defined as the number 
of 6-h track points per unit area, calculated by counting 
the track points in every 1° grid box and dividing by the 
area of the grid box. The track densities for the two 
ENSO states for each version of the model are shown in 
Fig. 8. The results show large similarities between the 
simulations from the Genesis-Only (Figs. 8a, d) and Both 
(Figs. 8c, f) model versions: the southeastern peak in 
density during El Nino moves to the northwest during 
La Nina. This shift in track density is consistent with the 
increased La Nina landfall rates on Vietnam, China, the 
Korean Peninsula, and Japan seen in Fig. 7c. The density 
shift is a consequence of shifting genesis regions, and is 
not observed in the Track-Only model (Figs. 8b,e). 
However, the Track-Only model shows higher track 
density overall in either ENSO state, with a notable 
northward increase in track density, especially during El 
Nino states. 
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GENESIS-ONLY 





DISTANCE ALONG COAST IN 1000KM 

Fig. 7. Ensemble-mean landfall rates for 1-yr model runs for ENSO indices of +2 (red) and 
-2 (blue) are shown for the coasts of Japan, the Philippines, and the two southernmost regions 
of the main coast. Each panel represents a different model configuration: (a) Genesis-Only 
(only ENSO dependence in the genesis component), (b) Track-Only (only ENSO de- 
pendence in the track-propagation component), and (c) Both (ENSO dependence in both), 
(d) The landfall rate difference between the two ENSO states from the Both model. The 
shading in (c) and (d) indicates the jackknife uncertainty at the 95% confidence level. The 
coastal segments with significant increases in La Nina year landfall are indicated in blue in (c). 


4. Discussion 

We have constructed a statistical track model for TCs 
in the WNP with seasonal climate and ENSO depen- 
dence in the genesis and track-propagation components 
with the goal of estimating the ENSO influence on East 
Asian landfall rates. The model implements local re- 
gression using seasonally varying climate variables and 
the JASO Nino-3.4 index as predictors. The advantage of 


a basinwide track model such as this compared with 
analysis of landfalls directly is that it brings to bear orders 
of magnitude more data to landfall estimation, essentially 
projecting information from the entire basin onto coastal 
regions of interest for climate states of interest. 

Stochastic simulations of TCs in extreme ENSO states 
reveal that ENSO affects landfall locations and fre- 
quency. We compute mean landfall rates and uncer- 
tainties about the means for strong El Nino and La Nina 
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Fig. 8. Track densities in units of 6-hourly track points per square kilometer for the three different versions of the model for (a)-(c) 
El Nino and (d)-(f) La Nina years. Note the effects of the genesis region shift for the (left) Genesis-Only and (right) Both versions 
and the high track density in the South China Sea for the (center) Track-Only version that agrees with the higher landfall rates in 
simulations. 


states. There are regions of significantly higher landfall 
rates in the Korean Peninsula, China, Vietnam, Japan, 
and the northern Philippines during La Nina. The majority 
of the ENSO-landfall relationship is mediated by the 
influence of ENSO on genesis rates and geographic 
distribution. However, ENSO influence on TC tracks 
adds and subtracts from the genesis-landfall effect in 
various coastal regions. 

The model is also evaluated by comparing simulated 
landfall rates to direct historical landfalls. Overall, the 
model performs well: in most regions historical landfall 
rates are inside the spread of model results across large 
ensemble simulations, indicating a lack of bias. There 
are a few regions where the model has biased landfall 
rates. These biases have been reduced, however, com- 
pared with the YH11 model version, which had no sea- 
sonality and no ENSO-track effects. The performance 
in historical simulations also indicates that the model 
has skill in simulating TC landfalls for all ENSO states, 
including neutral and weak events. The fact that our 
model has skill while using relatively few predictors 
highlights the importance of ENSO and seasonal effects 
on WNP TCs. 


Future directions for this work include developing 
a model component for the intensity and size of each TC, 
similar to Vickery and Twisdalc (1995) and Vickery et al. 
(2000), to make a more complete risk analysis. Studies 
such as Irish et al. (2008) and Weisberg and Zheng (2006) 
show that it is the complex relationship between fac- 
tors such as translational speed, storm size, intensity, 
bathymetry, and orientation that determines damage. 
Further, modeling of other TC-climate relationships may 
also be explored in the future where the model is dependent 
on other climate states. Advances in the understanding 
of the relationships between WNP TC activity and the 
Madden-Julian oscillation, quasi-biennial oscillation (Ho 
et al. 2009), central Pacific SST warming events (Kim et al. 
2011), and even the more remote Antarctic Oscillation 
(Ho et al. 2005) and Arctic Oscillation (Choi and Byun 
2010) motivate their additional consideration as statistical 
predictors to further improve the accuracy of our model. 
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